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We study how universal properties of quantum quenches across critical points are modified by a 
weak coupling to thermal dissipation, focusing on the paradigmatic case of the transverse field Ising 
model. Beyond the standard quench-induced Kibble-Zurek defect production in the absence of the 
bath, the bath contributes extra thermal defects. We show that spatial correlations in the noise 
produced by the bath can play a crucial role: one obtains quantitatively different scaling regimes 
depending on whether the correlation length of the noise is smaller or larger than the Kibble-Zurek 
length associated with the quench speed, and the thermal length set by temperature. For the case of 
spatially-correlated bath noise, additional thermal defect generation is restricted to a window that is 
both quantum critical and excluded from the non-equilibrium regime surrounding the critical point. 

We map the dissipative quench problem to a set of effectively independent dissipative Landau-Zener 
problems. Using this mapping along with both analytic and numerical calculations allows us to find 
the scaling of the excess defect density produced in the quench, and suggests a generic picture for 
such dissipative quenches. 

PACS numbers: 64.60.Ht, 05.30.Rt, 03.65.Yz 


Introduction- Quantum quenches involve the explicit 
time-dependent tuning of a Hamiltonian and are among 
the most basic and generic of phenomena in quantum 
many-body dynamics. They have garnered consider¬ 
able recent interest, particularly with regards to tuning 
through a quantum critical point, for which the quan¬ 
tum Kibble-Zurek (KZ) mechanism forms a standard 
paradigm M, analogous to the corresponding argu¬ 
ment for classical, thermal quenches [THlOj. The basic 
idea is that for gapped quantum phases separated by a 
gapless critical point, as one approaches the critical point, 
the system’s relaxation time (i.e. inverse gap, A“^) di¬ 
verges. When this relaxation time becomes comparable 
to the time associated with the quench speed v, the sys¬ 
tem falls out of equilibrium. This non-equilibrium regime 
dictates post-quench behavior, such as deviations from 
the final ground state and corresponding defect densities. 
One obtains universal scaling behavior of such quantities, 
ultimately governed by the proximity to the quantum 
critical point. 

In this work we turn our attention to the important 
yet relatively unstudied effect of thermal dissipation on 
the quantum Kibble-Zurek scenario. We focus on the 
prototypical system of a quench in a transverse-field 
Ising model (TFIM), where a dissipative thermal bath 
produces noise in the transverse field on each lattice 
site. Our study reveals that such dissipative quantum 
quenches are strongly influenced by the interplay between 
different fundamental length scales. 

A crucial role is played by the ratio of the bath-noise 
correlation length to both the length scale of the 
dissipation-free Kibble Zurek problem ^kz oc 1/\/u and 
to the thermal length scale oc l/k^T. Patane et 


al. [ill Hi investigated the limit ^ ^kz,Ct, where 
spatial noise correlations are essentially irrelevant, and 
each lattice site is effectively coupled to an independent 
bath. In contrast, we analyze the opposite but equally 
important limit where the noise correlation length, while 
finite, is nonetheless larger than both ^kz and In 
this case, the long-range bath noise correlations lead to 
strikingly different behaviours. 

For this correlated noise regime and generic form of 
system-bath coupling, our work suggests the following 
general picture, as graphically depicted in Fig.[T] During 
the quench, thermal effects come into play only under 
certain conditions. First, at a given time t, the bath 
temperature needs to be large enough to be able to pro¬ 
duce defects, implying that one must be in the quantum 
critical regime ksT > A{t) [^. However, this is not 
enough: even if the temperature is large enough, the co¬ 
herent system dynamics must be able to equilibrate in 
order for the weak dissipation to also play a role. This 
means that thermal defect generation is suppressed in the 
“non-equilibrium” Kibble-Zurek regime close to the crit¬ 
ical point. Thus, dissipation only gives rise to additional 
defects during (at most) a limited portion of the quench 
protocol. 


Focusing on temperatures small enough that only long- 
wavelength excitations can be produced, the above argu¬ 
ments lead to a general prediction that a weakly-coupled 
bath only gives rise to additional defect generation (above 
the standard zero-temperature Kibble-Zurek prediction) 
if the temperature satisfies ksT > ksTj^in- As detailed 
in the EPAPS 14|,we find that the Kibble-Zurek cri¬ 
terion for the cross-over scale into the non-equilibrium 
regime combined with constraints for being in 
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FIG. 1: (Color online) Quench trajectories for two temper¬ 
atures Ti and T 2 at a fixed quench rate. At Ti, far from 
the quantum critical point (QCP), the spectral gap A is too 
large for thermal excitations to be created. Thermal defect 
generation thus turns on only once the trajectory enters the 
quantum critical regime ksT > A. It is however suppressed 
again once the system enters the non-equilibrium regime near 
the QCP, where the intrinsic relaxational rate is slower than 
the quench rate. The thick yellow dashed line indicates the 
portion of the quench where thermal defects are produced. 
For low temperature (T 2 ), the entire quantum critical regime 
is subsumed within the non-equilibrium regime near the QCP, 
and hence there is no appreciable thermal defect generation. 


the quantum critical regime result in the following scaling 
behavior: 


AsT^inOC (1) 

Here v is the critical exponent describing the divergence 
of the system correlation length near the transition, and 
z is the corresponding dynamical critical exponent. For 
the specific case of the transverse field Ising model (where 
= z = I), we find that the excess defect density scales 
as (k^T)^/v as long as y/v k^T, but is strongly sup¬ 
pressed for lower temperatures. 

To rigorously establish the above picture, we utilize a 
major technical advantage of large : one can map the 
dissipative quench problem to an ensemble of dissipative 
Landau-Zener problems [IM 3 . We justify this mapping 
in what follows, and then make use of both analytic and 
numerical approaches to calculate the production of de¬ 
fects in the quench in the presence of a weak coupling to 
spatially-correlated dissipation. 

Model- We consider a one-dimensional TFIM subject 
to noise in the transverse field via a coupling to a thermal 
bosonic bath. We primarily focus on the limit where the 
bath noise is spatially uniform, like the average transverse 
field itself. The net Hamiltonian takes the form H = 
Hs + Hb + HSB, where 

Hs + Hsb = -JY. - (^ + ^) E (2) 

j 3 

represents the TFIM system and system-bath coupling 
Hamiltonian, Hb = is the bath Hamiltonian 


and X = [bu + h.c. j . Here, denote Pauli ma¬ 

trices for the spin at site J is the exchange coupling, h 
a Zeeman field in the z-direction, and we have set h = 1. 
Without loss of generality, we take J,h>0. The system 
is in an Ising ferromagnetically ordered phase when h < 
J, while for h > J the system is in a paramagnetic phase. 
The two phases are separated by a quantum critical point 
at he = J■ The bath is characterized by the spectral den¬ 
sity — Wi/) = , where 7 

is the dimensionless coupling strength and lOc is a cut-off 
frequency. We focus on the standard case of an Ohmic 
bath, where s = 1 . 

Employing the standard Jordan-Wigner transforma¬ 
tion [20| , we re-express the Hamiltonian in terms of spin¬ 
less fermions, Ci = (n,<. d|) a'l. Working in momentum 
space and setting the lattice constant a = 1 , we have: 


Hs + Hsb 


E c-k\{Hk,s+ Hk,SB) 

0<fc<7r 



with 'Hk,s 


A^ -a. 


Hk,SB 


X 0 

0 -w 


, ( 3 ) 


and where = 2 /i — 2 Jcos(fc) and Aj, = 2 Jsin(fc). 
The corresponding energy dispersion of Hs is given by 
Ck = yj^k + |Afcp. Critical exponents z/ = z = 1 
can be extracted from critical gap behavior, namely, 
Afc=o ^ \h — hc\ (near h = he = J) and Afc ~ |A:| (for 
/c — > 0 when h = he). 

Mapping to the Landau-Zener problem- We briefly re¬ 
call the standard treatment of quench dynamics in the 
absence of dissipation [l|. Each term in Hs (c.f. Eq. ©) 
describes an effective two-level system, where the two 
states correspond to having the orbitals {k, —k) either 
both empty or both occupied. The most common quench 
protocol involves a linear ramp of the form 


h{t) = —vt. 


( 4 ) 


For concreteness, we consider a quench that starts in a 
paramagnetic phase at t = —00 and ends in the ferro¬ 
magnetic phase at t = 0. With this time dependence, 
each term in Hs describes a Landau-Zener problem [il- 
[H: a two level system subject to a constant x mag¬ 
netic field (Afc) and linearly time-varying z magnetic 
field (^fc = —2vt — 2J cosk), which is taken through the 
avoided crossing occurring at ^k = 0 . 

Evaluating the density of defects ukz produced by the 
quench now amounts to calculating the final population 
of quasiparticles at the end of the quench. For slow 
quenches (with velocity u J^), the dominant contribu¬ 
tion to this population comes from low-fc fermion modes, 
where the gap Afc at the avoided crossing is the smallest. 
As discussed in Ref. 0 , the excitation probability of such 
modes is well-described by the asymptotic, infinite-time 
LZ probability that the two-level system (TLS) transi- 












3 


tions to the final excited state [BSil: 

Pk ~ exp {-ttAI/v) = exp {-irzk ), 


(5) 


where we have introduced the adiabaticity parameter 
Zk = A'^/v. Using this expression and integrating over 
all momentum modes, one obtains the power-law form 
nKZ consistent with the Kibble-Zurek form 

Mapping to dissipative Landau-Zener problem- Upon 
including the uniform coupling to the bath, we see from 
Eq. ([3|) that the translational invariance of the fermionic 
system is maintained. Thus, for each /c > 0, we have 
an effective two-level system whose detuning is linearly 
ramped in time and fluctuates now. The external noise 
X, however, couples to every k mode. The single-particle 
fermion Green functions for each momentum fc, however, 


are decoupled from one another (see [1^ for details). As 
a result, the defect density can be rigorously calculated 
by assuming that each k mode is coupled to its own, 
independent dissipative bath. This makes the mapping 
of our dissipative TFIM onto an ensemble of dissipative 
Landau Zener systems complete for the quantity of in¬ 
terest. Such dissipative LZ problems with diagonal noise 
have been well-studied in the literature 2ll429l|. 

Analytic estimates- The above mapping enables us to 
invoke results for the dissipative LZ problem, first estab¬ 
lished by Ao and Rammer 2l|, [22| , to evaluate the contri¬ 
bution of thermal dissipation towards defect production 
during the quench. For modes having adiabaticity pa¬ 
rameter Zk ^ I (c.f Eq. ([5|)), the transition through the 
avoided crossing is not adiabatic, and there is a large 
probability for the effective TLS to be excited even with¬ 
out dissipation. For such modes, dissipation (to lead¬ 
ing order) does not yield any additional probability of 
excitation beyond the LZ expression in Eq. Note 

that dissipation can only give rise to transitions between 
eigenstates which are coherent superpositions of (Tz eigen¬ 
states. While the instantaneous eigenstates have this 
form near the avoided crossing, for the fast modes with 
Zk ^ I, there is not enough time for these eigenstates to 
form. 

In contrast, for modes having Zfe > 1, the transition 
through the avoided crossing is very nearly adiabatic; 
without dissipation, they remain close to the ground 
state. For such “slow” modes, the system spends enough 
time near the avoided crossing for dissipation to give rise 
to transitions to the excited state. This leads to an addi¬ 
tional dissipation-induced excitation probability SPk that 
for weak dissipation simply adds to the LZ expression in 
Eq. ([5|). The probability SPk was explicitly calculated 
in Refs. 2l|, l22| and has a simple form consistent with a 
Golden rule calculation: 

SPk=‘^^Ji2Ak)nB[2Ak,T], ( 6 ) 

where nB[E,T] is a Bose-Einstein distribution evaluated 
at energy E and at the temperature T of the dissipative 


bath. This is simply the Fermi’s Golden Rule rate for 
excitation of the TLS at the avoided crossing via absorp¬ 
tion of a bath phonon (Texc = 2TTj{2Ak)nB{2Ak)), mul¬ 
tiplied by the effective time spent at the avoided crossing 

(Aross ~ Ak/v). 

Hence, as sketched in Fig. [U the dissipation only in¬ 
fluences modes where Ak is sufficiently large to yield a 
near-adiabatic transition, but also small enough that the 
Bose-Einstein factor in Eq. m is appreciable (i.e. small 
enough that the bath temperature is sufficient to gen¬ 
erate an excitation). Thus, the only modes affected by 
dissipation satisfy: 

{Akz = Vv) ^ Ak < ^bT. (7) 

For sufficiently slow quench velocities, all the relevant 
modes satisfying this condition correspond to small 
k, where Ak ^ Jk. This justifies the general pic¬ 
ture provided in Fig. [TJ Outside the quantum criti¬ 
cal regime, A = > T, and there are no modes 

which satisfy Eq. 0 ; dissipation thus has no effect here 
(i.e. blue/purple region in Fig. [I}. In the non-equilibrium 
regime (orange in Fig. [1]), A^^q < and v is too 

large for dissipation to affect any of the inodes. Fi¬ 
nally, in the portion of the quantum-critical region that 
lies outside the non-equilibrium regime, we have instead 
Akz < Afe=o < ^bT; we thus have modes satisfying 
Eq. ([7|) and dissipation can create excitations. Here, we 
observe dynamics different from the KZ scenario. 

Integrating SPk of Eq.® over the range defined by 
Eq. ([7]), and assuming k^T fyu, we find that the total 
dissipation-induced excitation density obeys: 

nth X^kBT[{kBTf/v]. (8) 


This result is consistent with a more general ansatz for 
a scaling form for the dissipation-induced defect density, 
see EPAPS Q. 

Numerics - To confirm the picture presented above, we 
perform a rigorous numerical analysis of a TFIM quench 
within the dissipative LZ framework. We employ weak 
system-bath coupling techniques where simple but ap¬ 
proximate Markovian approaches, specifically nonequi¬ 
librium Bloc h eq uations (NBEs) [30|, have been shown 
to be reliable 13 III in comparison to numerically-exact ap¬ 
proaches [^, l27l | (details in 0 ). 

The quench is explicitly started at hit^) = — 8 J and 
ended at h{0) = 0, implying v = —8J/to. We have 
ensured that results do not depend on this initializa¬ 
tion. Solving the adiabatic-Markovian master equation 
for each k, we obtain the total probability P{k^v;"f,T) 
for a given effective TLS to end in the excited state. The 
thermal contribution Pth(fc, v; 7 , T) is determined by sub¬ 
tracting the dissipation-free probability P{k, u; 7 = 0, T) 
from the total one. Integrating over all k we find the 
total excitation density produced during the quench. 


ntot('^;7^^) = dk P{k,v;j,T) = nKz{v)-\-nthiv;j,T), 
Jo 
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FIG. 2: (Color online) Total excitation density ntot(n; 7 , T) 
versus quench speed v for various temperatures, ujc J, and 
a system-bath coupling 7 = 1.3-10“®. The red full line shows 
nKz(v), i.e. the excitation density for 7 = 0 . 


where tikz is the defect density without dissipation, and 
nth is the additional defect density due to the thermal 
dissipation. 

In Fig. M we plot the total defect density produced by 
the quench as a function of quench velocity v; different 
curves correspond to different bath temperatures. The 
defect density for the dissipation-free case 7 = 0 is also 
plotted (red solid line); this curve exhibits the standard 
y/v Kibble Zurek scaling. For finite 7 and T we observe 
with decreasing v an increasing defect probability due to 
additional thermal defects. Thus, at a fixed temperature 
a minimal total defect density is observed for an optimal 
quench speed. The thermal defect density also increases 
with increasing temperature. 


To focus on the bath contribution to the defect produc¬ 
tion during the quench, in Fig. |3]we now plot results for 
the excess thermal defect density as a function of scaled 
temperature. The analytic estimate predicts a scaling 
nt\i {UbTY/ v for k^T ^ J (see Eq. ([U) which 

is indeed observed (see data and orange sold line in Fig. 
El). For smaller T, we see that the thermal defect den¬ 
sity is suppressed, also in accordance with our prediction 
that thermal defect production is suppressed in the “non¬ 
equilibrium” regime 1/A(f) > \t\. 

The inset of Fig. El plots the momentum-resolved ther¬ 
mal defect density Pth(fc, u; 7 , T) exhibiting a clear peak. 
With increasing temperature the momentum of the peak 
maximum increases in line with expectation from the be¬ 
haviour of the minimum in the adiabatic regime of the 
dissipative Landau-Zener dynamics 26|, l27|. For temper¬ 
atures T —>■ y/v the main contribution is from modes at 
k ~ y/v. Note that here the adiabatic - Markovian mas¬ 
ter equation is least reliable 3l| and thus the observed 
behaviour in the low temperature regime in Fig.|31cannot 



FIG. 3: (Color online) Excess thermal defect density associ¬ 
ated with the quench, as a function of scaled temperature, 
for three different quench velocities v. For ^/v ^ T < J, 
the data follows nth <x T^/v in line with the scaling form 
given in Eq. (l 8 l (orange sold line). Inset: Plot of momentum- 
resolved thermal defect density Pth{k,v,y,T) for quench ve¬ 
locity V = 1.4 • 10“® and three temperatures. 


be taken to be conclusive. 

Finite bath correlation length- While the above anal¬ 
ysis was for noise that is correlated over the entire lat¬ 
tice, it remains valid for situations where the bath noise 
has a finite spatial correlation length as long as 
this length is much longer than all other relevant length 
scales. In particular, this requires be larger than both 
the Kibble-Zurek length ^kz = Jafy/v and the thermal 
length = Ja/ksT. It does not however require that 
approach the size of the system. The fact that the 
system correlation length diverges at the QCP is not rel¬ 
evant; heuristically, the bath only plays a role outside the 
non-equilibrium regime surrounding the QCP, and thus 
does not see the full divergence of this length. A more 
rigorous justification for why our results apply whenever 
Cn ^ CkZjCt (based on a Keldysh analysis) is provided 
in the EPAPS [ 3 . In the opposite limit of a very small 
correlation length “C ■Cr, one approaches the limit of 
spatially-uncorrelated dissipation, as studied by Patane 
et al [ll|, [I2| , with a very different temperature depen¬ 
dence of the thermal defect density from Eq. We thus 
see that the ratios /■Ckz and / £,t will in general play 

a crucial role in determining the influence of dissipation 
on a quench. 

Conclusions- We have studied the general problem 
of how thermal dissipation and noise give rise to defect 
generation in a quantum quench beyond the amount pre¬ 
dicted for zero dissipation by the Kibble-Zurek scenario. 
We have explicitly analyzed quenches in the TEIM cou¬ 
pled to a global thermal bath by mapping to a set of 
effectively independent dissipative Landau-Zener prob¬ 
lems. By analytical estimates as well as numerical cal¬ 
culations, we have shown that the excess thermal defect 
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density scales as {k^TY/v as long as -^k^T J in 
clear contrast to previous results obtained for spatially- 
uncorrelated dissipation. These first steps toward under¬ 
standing the interplay between quench and dissipation 
elicit a host of further studies, for instance, regarding 
more complex systems, different crossover regimes, gen¬ 
eral scaling laws, and effect of alternative bath or cou¬ 
pling forms. As with non-dissipative quantum quenches, 
certain cold atomic and magnetic systems would form the 
experimental testbeds for these studies. 

We thank Alessandro Silva and Anatoli Polovnikov for 
useful conversations. PN acknowledges financial support 
by the excellence cluster “The Hamburg Centre for Ultra¬ 
fast Imaging - Structure, Dynamics and Control of Mat¬ 
ter at the Atomic Scale” and by the project “Nonequi¬ 
librium Quantum Energy Transport through Nanostruc¬ 
tures” (NA 394/2-1) of the Deutsche Forschungsgemein- 
schaft. SV acknowledges the National Science Founda¬ 
tion grant DMR 0644022-CAR. AAC thanks NSERC. 
AAC and SV also thank the Kavli Institute for Theo¬ 
retical Physics for their hospitality. 


[1] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- 
tore, Rev. Mod. Phys. 83, 863 (2011). 

[2] J. Dziarmaga, Phys Rev Lett 95, 245701 (2005). 

[3] A. Polkovnikov, Phys. Rev. B 72, 161201 (2005). 

[4] W. H. Zurek, U. Dorner, and P. Zoller, Phys. Rev. Lett. 
95, 105701 (2005). 

[5] J. Dziarmaga, Advances in Physics 59, 1063 (2010). 

[6] A. Dutta, U. Divakaran, D. Sen, B. K. Chakrabarti, 
T. F. Rosenbaum, and G. Aeppli, ArXiv e-prints (2010), 
1012.0653. 

[7] T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 

[8] T. W. B. Kibble, Phys. Rep. 67, 183 (1980). 

[9] W. H. Zurek, Nature 317, 505 (1985). 

[10] W. H. Zurek, Phys. Rep. 276, 177 (1996). 

[11] D. Patane, A. Silva, L. Amico, R. Fazio, and G. E. San¬ 
toro, Phys Rev Lett 101, 175701 (2008). 

[12] D. Patane, L. Amico, A. Silva, R. Fazio, and G. Santoro, 
Phys. Rev. B 80, 024302 (2009). 

[13] S. Sachdev, Quantum Phase Transitions (Gambridge 
University Press, 2011), 2nd ed. 

[14] See supplementary information for further details. 
(2015). 

[15] A. Ghandran, A. Erez, S. Gubser, and S. L. Sondhi, Phys. 
Rev. B 86, 064304 (2012). 

[16] L. D. Landau, Phys. Z. Sowjetunion 1, 89 (1932). 

[17] G. Zener, Proc. Roy. Soc. London A 137, 696 (1932). 

[18] E. G. C. Stueckelberg, Helv. Phys. Acta 5, 369 (1932). 

[19] E. Majorana, Nuovo Cimento 9, 43 (1932). 

[20] P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

[21] P. Ao and J. Rammer, Phys. Rev. Lett. 62, 3004 (1989). 

[22] P. Ao and J. Rammer, Phys. Rev. B 43, 5397 (1991). 

[23] E. Shimshoni and Y. Gefen, Ann. Phys.-New York 210, 
16 (1991). 

[24] Y. Kayanuma and H. Nakayama, Phys. Rev. B 57, 13099 
(1998). 


[25] M. Wubs, K. Saito, S. Kohler, P. Hanggi, and 
Y. Kayanuma, Phys. Rev. Lett. 97, 200404 (2006). 

[26] P. Nalbach and M. Thorwart, Phys. Rev. Lett. 103, 
220401 (2009). 

[27] P. Nalbach and M. Thorwart, Ghem. Phys. 375, 234 

( 2010 ). 

[28] P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 
82, 032118 (2010). 

[29] P. P. Orth, A. Imambekov, and K. Le Hur, 

Phys. Rev. B 87, 014305 (2013), URL 

http://link.aps.org/doi/10.1103/PhysRevB.87.014305 

[30] P. Nalbach, J. Knorzer, and S. Ludwig, Phys. Rev. B 87, 
165425 (2013). 

[31] P. Nalbach, Phys. Rev. A 90, 042112 (2014). 

[32] A. Kamenev, Field Theory of Non-Equilibrium Systems 
(Gambridge University Press, 2011). 

[33] W. H. Press, S. A. Teukolsky, W. T. 

Vetterling, and B. P. Flannery, 

Numerical Recipes: The Art of Scientific Computing 
(Cambridge University Press, 2007), 3rd ed. 


APPENDIX - SCALING ARGUMENTS 

Here we generalize the scaling arguments presented in 
previous works Qim with regards to the universal 
non-equilibrium dynamcs associated with Kibble-Zurek- 
type physics. We focus on the influence of dissipation 
from a weakly coupled thermal bath which produces 
spatially-correlated noise. We also focus on a system- 
bath coupling which is analogous to transverse-field noise 
in the transverse-field Ising model, i.e. the system bath 
coupling Hamiltonian commutes with the coherent sys¬ 
tem Hamiltonian away from the critical point. 

As the quench process evolves through the quantum 
critical point, we expect excitations produced by cou¬ 
pling to the thermal bath, nth-, to become important 
in the vicinity of the critical point where the system’s 
energy gap A is small compared to the bath tempera¬ 
ture scale ksT. However, very close to the critical point, 
the system’s relaxational time scale, r, diverges and the 
quench process becomes too rapid for the system to form 
the coherences needed to allow bath-induced transitions. 
These constraints provide the following bounds for the 
regime in which thermal excitations are significant and 
the following scaling forms for their behavior in terms of 
characteristic energy/time scales. 

Consider a quantum phase transition characterized by 
a parameter a such that the quantum critical point oc¬ 
curs at a value Oc [l^. Close to the critical point, the 
system’s typical correlation length scale diverges as 

e (9) 

where v is the associated critical exponent and <5 = ja — 
ad measures the deviation from the critical point. The 
relaxational time diverges as 

r ^ 


g-i.z 


( 10 ) 
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where z is the dynamic critical exponent. 

With regards to Kibble-Zurek scaling behavior[l[ 
[l5l| . consider a linear quench at a characteristic rate v~^, 

a{t) = ttc + vt. (11) 

A given deviation 5 thus occurs at a time t{5) = S/v; 
|t((5)| represents the time remaining until the system 
reaches the critical point. Now the quench enters the 
non-equilibrium regime under the condition 

t{6) < T, (12) 

or, from Eq. (flUl) . 5/v < 5~''^. This cross-over criterion 
establishes a scaling relationship between the quench rate 
and distance to criticality, namely 

5^+''^ ~ V. (13) 


From this relationship, one can derive the scaling of stan¬ 
dard Kibble-Zurek variables. For instance, the density 
of defects produced in the non-equilibrium region scales 

as ~ ^ yVd/{l+vz) ^ gj: ^Jj^g ggg_ 

dimensional transverse Ising system, where v = z = d = 
1 , one obtains the well known scaling behavior no ^ y/v. 

With regards to the regime of interest here, namely, 
that of thermal excitations, the two constraints men¬ 
tioned above need to be satisfied. First, as discussed 
previously in Ref. 0,[il , the bath temperature must be 
greater than the gap, i.e. fc^T > A, yielding the scaling 
relationship ksT ~ 5''^. Second, the quench must fall 
short of entering the non-equilibrium regime, or equiva¬ 
lently, t{5) > {kBT)~^. This condition, combined with 
the scaling relationship obtained from the first constraint 
and with Eq. 0 provides the following temperature 
lower bound for a given quench rate: 

> u, (14) 


as presented in Eq. (1) in our main text. 

In Ref. 1^, similar to equilibrium quantum critical 
scaling, the effect of finite temperature has been dis¬ 
cussed in the context of non-equilibrium quantum critical 
scaling. This is done by introducing a dimensionless pa¬ 
rameter that is naturally defined by the scaling relation¬ 
ship between temperature and quench rate in Eq. 0. 
It is interesting to note that the excess thermal density 
that we predict in Eq. (8) of the main text, given that it 
has dimensions of inverse volume, is consistent with the 
scaling form hypothesized in Ref. for generic situa¬ 
tions: 


nth - (15) 

where is a scaling function. Specifically, in the one¬ 
dimensional transverse Ising case, our arguments show 
that thermal excitations are important in the regime 
(fcsT’)^ > V and that they respect the form nth 
kBTF[kBT/^. 


We emphasize however that care needs to be taken in 
appl ying the above scaling. The situation presented by 
Ref. [1^ presents a closed system having an initial tem¬ 
perature fc^T, and in general, unlike our case, thermal 
effects need not give rise to separate additional contribu¬ 
tions above the zero-temperature Kibble-Zurek contribu¬ 
tions. We believe that adherence to the expected form 
is tied to the weak nature of the bath coupling, as well 
as the particular choice of an Ohmic bath spectral func¬ 
tion. Considering the effects of more general bath spec¬ 
tral functions and stronger couplings would make for an 
interesting and challenging study. 


APPENDIX - MAPPING TO INDEPENDENT 
DISSIPATIVE LANDAU-ZENER TRANSITIONS 


In the main text, we argue that in the limiting case 
where each site of our TFIM couples to the same dissipa¬ 
tive bath (a “global” system-bath coupling), we can map 
our dissipative quench problem onto a set of independent 
dissipative Landau-Zener problems. One might worry 
that this mapping is only approximate, as it ignores cor¬ 
relations between different fermionic modes induced by 
the globally-coupled bath. While such correlations will 
play a role for some physical observables, they play no 
role in determining single-particle properties, such as the 
quasiparticle occupancies which we focus on. This is a di¬ 
rect consequence of the fact that the system-bath Hamil¬ 
tonian of Eq. (3) in the main text conserves the momen¬ 
tum of the fermionic system. 

More formally, the occupancy of a given quasparticle 
mode with momentum k can be written 

Pk{t) = (4Wcfc(i)) = (16) 


where we have introduced the standard lesser Keldysh 
Green’s function associated with this mode [s^- If one 
treats the coupling to the bath as a perturbation of 
the coherent (dissipation-free) system, the lesser Green’s 
function appearing above is completely determined by 
the Keldysh self energies I]^(t,f') associated with the 
system-bath interaction; here, the index a can take the 
values i?. A, K, corres pon ding to retarded, advanced and 
Keldysh self-energies 32|. Note the self-energy must be 
diagonal in momentum, as electronic momentum is con¬ 
served. Gonsider an arbitrary self energy diagram for 
As fermion momentum is conserved at each 
system-bath interaction vertex, all internal fermion prop¬ 
agators in this diagram involve the same momentum k. 
Heuristically, this implies that at least for single-particle 
Green’s functions, a given fermion mode with momentum 
k does not know about other modes having a different 
momentum k' ^ k. 

It follows that we would obtain exactly the same dia¬ 
grammatic expansion (and hence result for (t, t)) if we 
had coupled each fermion mode to its own independent 
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bath. Formally, this means modifying the system-bath 
Hamiltonian in Eq. (3) of the main text as follows: 

Hsb c-fe] 

0<fc<7r 

(17) 

We now have an independent bath for each k mode, with 
a corresponding noise operator Xk = ^k,v{bk,v+h.c.). 

Each of these baths (labelled by k) has identical proper¬ 
ties to the bath appearing in our starting Hamiltonian. 
They all have the same temperature T and identical spec¬ 
tral densities: Jk{^) = where 77(w) is the spectral 

density of our original bath (noise operator X), as given 
after Eq. (2) in the main text. 

Thus, for computing quasiparticle occupancies, we can 
exactly treat each fermion k mode as being effectively 
coupled to its own independent bath. This then rigor¬ 
ously justifies our mapping to an ensemble of uncoupled 
dissipative Landau-Zener problems. 


Xk 0 
0 -Xk 









APPENDIX - ADIABATIC MARKOVIAN 
MASTER EQUATION 

Letting \j{t)) {j = 1,2) denote the eigenstates 
of the instantaneous coherent Hamiltonian 'Hk,s{t) in 
Eq. (3) (in the main text) (with eigenergies ztEk{t) = 

±\J{ik{t))^ + ), the matrix elements of the ef¬ 
fective statistical operator pk(t) of the Landau-Zener 
system after tracing out the bath degrees of free¬ 
dom are parametrized in terms of a 3-vector rk{t) as 
(j(0l/3fe(^)l/(0) = - U- • ^)h': where r is the vector 

of Pauli matrices. Suppressing the k index for clarity, 
the adiabatic Markovian master equation takes the form 

HI: 

/-7i(t) 0 0(t) \ 

dtr=\ 0 - 72(0 -2E{t) r-f- 7 i(t)req, (18) 

V -m 2E{t) -72(t) ) 

where 9{t) = arctan[^fc(t)/Afc] and Feq = 

(tanh[/3£'(t)], 0,0). The terms proportional to 0 
describe coherent non-adiabatic evolution, while the 
time-dependent relaxation and dephasing rates are 

71 (t) = cos^ {9{t)) S[2E{t)] (19) 

72(i) = 57i(i) +sin^(fi'(i))'5[0] (20) 

where .S[u;] = 27r^7[a;] coth/3w/2 is the symmetrized spec¬ 
tral density of the bath noise. 

For a given quench protocol, i.e. ^k(t), we solve the 
adiabatic-Markovian master equation ([1^ using a stan¬ 
dard fourth order Runge Kutta scheme (33[ . Thus we get 
the time-dependent statistical operator for each k and, 
in turn, the probability P{k, v; 7 , T) for a given effective 


TLS to end up in the excited state at the end of the 
quench. To obtain the total excitation density produced 
during the quench, we integrate over all k using a stan¬ 
dard Romberg scheme [ 3 ^ . 


APPENDIX - EINITE BATH SPATIAL 
CORRELATIONS 


As discussed in the main text, our results for a 
spatially-uniform dissipative bath remain valid in the 
case where the bath noise has a finite spatial correlation 
length ,^ 7 v, as long as this length is much larger than the 
Kibble-Zurek length ^kz; one does not need to ap¬ 
proach the size of the system. To make this precise, we 
generalize the system-bath coupling so that there is a dis¬ 
tinct bath noise operator on each site, Hsb = J2j ^j^j- 
As usual, we take the bath to be an infinite collection of 
harmonic oscillators in thermal equilibrium, and take the 
Xj to be linear in the bath creation and destruction op¬ 
erators. We also take the bath to be in a translationally 
invariant state (unlike the model presented in Ref. flih . 
We consider a generic situation where (like the main text) 
the on-site noise is still described by an Ohmic spectral 
density but where the noise correlation decays ex¬ 

ponentially with distance 


{{X,[u;],Xk[co']}) 
2 'k5{uj -b w') 


= J[oj] coth 


w 


2 fcsT 


o-b-fcl/Cw 


( 21 ) 

The above form implies that we can express the 
Fourier-transformed bath noise operators as: 


E = A, A. (6,,. + 6%,,) (22) 

^ l V 


where the hq^v describe Einstein phonons with en¬ 
ergy uji,[q\ = Wi/. We use a normalization such 
(l/V)^^|Aqp = 1. The spectral function J7’[w] asso¬ 
ciated with the noise on any given site is then identical 
to that used in the main text. 

We next specialize to exponentially decaying spatial 
correlations, with a correlation length much smaller than 
the system size. This implies: 


|A,|^ 


2^n/ a 


(23) 


Unlike the global-bath model in the main text, the fi¬ 
nite correlation length here means that the bath can ex¬ 
change non-zero momentum with the system. We follow 
Ref.ll^ and use the weakness of the system-bath coupling 
to treat the system perturbatively, using a self-consistent 
Born approximation for the Keldysh self-energy of the 
fermion Green functions. Within this approximation, the 
Keldysh self energy for a fermion with momentum k is 












given self-consistently by 


Sfc(t2,172;, ti, (Ti) = (24) 

^ Gk+q{t2i <^2', ,ti,ai)D°{t2,<T2', ,ti,ai), 

q 

where CTj = ± denotes the forward and backward Keldysh 
contours, Gk [t, cr; t'a'] is a dressed fermion Keldysh Green 
function, and D^[t,(T\t'a'] is the unperturbed bosonic 
(equilibrium) Keldysh Green function for the bath op¬ 
erator Xq/Kq. In our model (where we assume Einstein 
phonons, corresponding to a frequency-independent ^n), 
this Green function is independent of q. 

We are now in a position to make estimates concern¬ 
ing the role of based on the behaviour of the imagi¬ 
nary part of the self-energies (which control bath-induced 
transitions). We will focus on transitions which are ther¬ 
mally enhanced, i.e. which involve the absorption or emis¬ 
sion of bath phonons having w < ksT. For a bath- 
induced scattering event taking a quasiparticle from mo¬ 
mentum k to k q, energy conservation and the fermion 
dispersion relation will determine the energy of the bath 
phonon involved. Assuming ksT <C J as always, this 
then naturally leads to the thermal length = aJ/ksT: 
the only transitions that are thermally enhanced involve 
momentum transfers with jgl < 1 /Ct- 

The simplest regime is where “C ^t- In this 
case, the only thermally-enhanced transitions have q <C 
1/^T 1/^Af, and the g-dependence of the structure fac¬ 

tor Aq plays no role: we can safely replace A, by Aq=o- 
In this case, the q integral in the self energy of Eq. (IMl) 
can be estimated as: 

n n /-i/Ct 

- \Aq\^Gk+q^-\Ao\^ Gk+q (25) 

d-Tr/a 27r J-l/^T 

cx {^N/a)T (26) 


In this limit where Aq can be treated as a constant, we 
recover the local-bath model studied in Ref. Il2, where 
there are no spatial correlations between bath noise op¬ 
erators Xj. This mapping to uncorrelated noise is valid 
irrespective of the value of ^n/'Ckz, where ^kz = aJ/y/v 
is the Kibble-Zurek length introduced in the main text. 
Note the explicit factor of T that emerges from the mo¬ 
mentum summation. 

Consider next the opposite regime, where ^ ^t- In 
this case, the structure factor Aq will suppress the con¬ 


tribution of large momentum transfers in the self-energy, 
as opposed to the bath temperature: the largest con¬ 
tributing |g| will be ~ I/Cat- H in addition we have 
^ ^KZi then we can also ignore the q dependence 
of the fermion propagator Gk+q in Eq. (IMl) . To under¬ 
stand this point, note that without dissipation, quasi¬ 
particle modes with k <C 1/Ckz will evolve diabatically 
during the quench (and become excited), while modes 
with k ^ 1/Ckz will evolve adiabatically (and hence re¬ 
main unpopulated). Correspondingly, for small fc, the 
quasiparticle modes k and k + q behave almost identi¬ 
cally when |g| < 1/Ckz, implying Gk ~ Gk+q- 

Thus, for ^ the q summation in the self¬ 

energy of Eq. (|M1) can be estimated as 



(27) 

(28) 


This is identical to having taken the global-coupling, 
^AT —>■ oo limit from the outset, i.e. having used 


Aq = ^S{q) 


(29) 


Thus, when is the largest length scale in the problem 
(i.e. ^Af ^ CkZi^t), the system does not know about the 
finite bath spatial correlation length and one gets 
the same results as for a model where —>■ oo. Note 

that in this large limit, the summation over q gives a 
temperature-independent result, in contrast to the small 
estimate in Eq. (l26l) . This difference is at the heart 
of why our global-coupling result for the thermal defect 
density scales as a lower power of temperature than the 
corresponding result found in Refs. luJ, [ij for a locally- 
coupled bath. 

Finally, the is the remaining case "C C Ckz- 
In this case, the bath correlation length cuts off large 
momentum transfers (as opposed to temperature). How¬ 
ever, this cutoff is not large enough to prevent coupling 
between very different fermionic modes, i.e. adiabatic and 
non-adiabatic modes. In this cross-over regime, neither 
the local model studied in Ref. [ll] nor global bath model 
studied in the main text are appropriate. 






